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Abstract 

We address the problem of antiferromagnetism in a two dimensional model of doped 
spin-Peierls system, at the classical and quantum levels. A Bethe-Peierls solution is 
derived for the classical model, with an ordering temperature proportional to the doping 
concentration. The quantum model is treated in a cluster renormalization group showing 
a finite randomness behavior and an antiferromagnetic susceptibility at low temperature. 
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1 Introduction 



The spin-Peierls transition at Tgp ~ 14 K in the inorganic quasi one dimensional spin-Peierls com- 
pound CuGeOs has attracted much interest fij. This transition is characterized by the appearance 
of a finite dimerization in the CuC>2 chains, and the opening of a spin gap. With CuGe03, it became 
possible to experiment the effect of doping in a spin-Peierls system. An antiferromagnetic (AF) 
phase was discovered upon replacing a fraction of the Cu ions (with S = 1/2) by magnetic ions with 
a different spin: Ni M (with S = 1) or Co || (with S = 3/2), or non magnetic ions: Zn [Q, [|, |6| 
or Mg Q. Also, the Ge sites (outside the Cu02 chains) can be substituted with Si leading to 
antiferromagnetism at low temperature. 

Recent experiments by Manabe et al. have shown the existence of a finite Neel temperature 
~ 25 mK, with a doping concentration as low as 0.12 % ||. The doping dependence of the Neel 
temperature obtained in these experiments suggests the absence of a critical concentration for the 
appearance of antiferromagnetism: at low doping the ordering temperature scales like InT/v oc 1/x Q. 

11 , ^] have focussed on the identification of the relevant low energy 



Early theoretical works 



degrees of freedom. Fukuyama, Tanimoto and Saito [11| have shown the coexistence between an- 
tiferromagnetism and dimerization in a doped spin-Peierls model. The degrees of freedom relevant 



to the low energy physics are solitonic spin-1/2 excitations pinned at the impurities [12]. These 
excitations are the building blocks of the theory in Refs. |l3| , 14, |l^]. These spin-1/2 objects interact 
via an exchange decaying exponentially with distance. Interchain interactions can be incorporated 
by considering the existence of a transverse correlation length, approximately one tenth of the lon- 
gitudinal correlation length, as recently proposed independently by Dobry et al. ]l6|, and Fabrizio, 
Melin and Souletie [15|. Numerical calculations with realistic spin-phonon couplings have provided 
a link between the microscopic Hamiltonian and the effective model of interacting spin-1/2 mo- 
ments [TO, |17], 19, [20]. The approach followed in Refs. |l^, 14, [l5| and continued in the present 
work relies on the treatment of disorder in the effective Hamiltonian. This allows to discuss the 
qualitative physics of large scale systems at a finite temperature, while the numerical methods have 
so far been limited to the ground state properties [p], 1C, 17, 18, 19, p0[. The scope of this article is to 



analyze the model beyond the percolation approximation used in Ref. 15], both at the classical and 
quantum levels. We first show in section || that the physics of the quantum Hamiltonian is already 
present in the classical Ising Hamiltonian, and give a rigorous derivation of mean field theory via a 
Bethe-Peierls treatment. The second purpose of the article is to show that the quantum Hamiltonian 
has an antiferromagnetic behavior at low temperature. The quantum model is treated in a cluster 
renormalization group (RG) calculation in sections || and || 
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2 The model 



We recall the model proposed in Ref. [13|. When impurities are introduced in a dimerized back- 
ground (for instance non magnetic impurities such as Zn), spin- 1/2 solitonic moments are released 
out of the dimerized pattern. These magnetic moments are pinned at the impurities due to interchain 
interactions [12]. This picture is in agreement with susceptibility experiments ||, 21, 22], indicating 
the release of one spin- 1/2 moment per Zn impurity at low doping. The interaction between two 
magnetic moments at a distance d originates from virtual excitations of the gaped dimerized back- 
ground, and decays exponentially with distance, with a characteristic length set by the correlation 
lengths £ x ~ 9 c along the chain direction (c-axis) p3| , 24], and £ y ~ £r/10 in the b-axis direction. 



These exchange interactions as well as the relevance of disorder were identified in Ref. [13] to play 
a crucial role in the establishment of antiferromagnetism. The low energy physics of a doped spin- 
Peierls system is represented by spin- 1/2 solitonic magnetic moments distributed at random with a 



concentration x on a square lattice, and interacting via a Heisenberg Hamiltonian [15, |R 

H 



£ J, 



(1) 



the exchange in Eq. |T] being staggered and decaying exponentially with distance: 

(_l)^+^+i Aexp 



J, 



d 



+ 



(2) 



0.1 x £ x the correlation length along 



with £ x ~ 9 c the correlation length along the c-axis and £ y 
the b-axis. Correlations along the a-axis are neglected. 

3 Classical Ising model 



We consider a model with Ising degrees of freedom, distributed randomly and interacting via the 
exchange in Eq. ||. The classical antiferromagnet has the same transition temperature as the classical 
ferromagnet. We consider therefore the ferromagnetic model to calculate the ordering temperature. 

3.1 One dimensional model 

A high temperature expansion leads to the exact form of the correlations in terms of a product over 
the bonds between the spins at sites and L: (<7o0"l) = Yl tanh (/JJ^j+i). We calculated numerically 
the disorder average to obtain the correlation length at a finite temperature. As shown on Fig. |l|, 
the average correlation length of the disordered model is larger than the typical correlation length 
£ = — l/[xln (tanh (/3T*))], with T* = Aexp (— l/(x£)) the exchange of the particular disorder 



realization where the magnetic moments are equally spaced. We calculated in Ref. [14] the correlation 
length of the quantum chain and found a similar result: the enhancement of the magnetic correlations 
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Figure 1: The average correlation length of the disordered Ising spin chain with nearest neighbor 
couplings is larger than the correlation length of an ordered system with the same concentration 
x = 0.01 above T*. We used A = 44.7 K, £ = 9. 

above T* due to disorder does not rely on the quantum nature of the coupling Hamiltonian in spite 
of a random singlet physics in the quantum chain, not present in the Ising chain. 

3.2 Determination of the exchange distribution 

We consider the exchanges to be drawn independently in a distribution P(J) resulting from the 
combination of randomness in the spatial distribution of the magnetic moments and exponentially 
decaying interactions, Eq. ||. The relevant exchanges are set by the spins the closest to each other. 
Therefore, given a spin at site (xq, yo), we need to determine the probability that one spin is found on 
the periphery of the ellipse [(x — xq)/^ x ] 2 + \{y — yo)/£,y] 2 = 7 2 > with no other spin inside the ellipse, 
and therefore an exchange Aexp (—7). We consider a system of total area A containing n spins. The 
probability to find no spin inside a subsystem of area 5 A is Pq = ^1 — ^\ ~ exp(— xSA), with 
x = n/A the doping concentration. Now the spacing distribution is 

P( 7 )=xL( 7 )exp(-x<5A(7)), (3) 

with L{pi) = d[A(^)]/d , y. In the one dimensional model, we have SA(j) = 2j£ x , and £(7) = 2^ x . In 
the two dimensional isotropic model with £ x = £ y = £, we have ^(7) = irj 2 ^ 2 , and £( 7 ) = 27r7£. 
In the quasi one dimensional model, SAfr) = 7r"f 2 £, x £, y , and L(-y) = 2iry£ a .l*y. The distribution P{j) 
of the isotropic and anisotropic two dimensional models is a Wigner distribution with a short scale 
"distance repulsion" . This repulsion will be shown not to affect the ordering properties. 
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3.3 Bethe-Peierls transition in the infinite coordination limit 



We now consider the Bethe-Peierls solution of the Ising model 25]. The lattice has a tree topol- 
ogy, with a forward branching ratio z — 1, and we calculate the magnetization of the site with the 
highest hierarchical level (top spin), in the presence of the other sites. We consider z — 1 trees and 
connect them to obtain a tree with one more generation (see Fig. ^). The recursion of the average 



magnetization of the top spin reads [26] 



x = n'=i(l + Y tanh Q3JQ) - Pj^l - Y t tanh {(3 J,)) 
n-=i(l + Y tanh (PA)) + n-=i(l - Y tanh {(3Ji)) ' 

where Yi, i = l,...,z — 1 are the magnetizations of the top spins with n generations, and X the 
magnetization of the top spin with n + 1 generations. We first consider the artificial situation where 
the ordering temperature is large compared to the exchange: Tb p 3> A, which turns out to be 
equivalent to assuming a large coordination. Eq. ||] is linearized into X = J2i=i Yi tanh (flJi), leading 
to the recursion of the magnetization ((X)) n+ i = (z — l)((tanh (j3J))){(X}) n , with the subscript 
n labeling the number of generations. This leads to the ordering temperature Tb p = (z — 1)((J}}, 
far above A if z 3> 1, and consistent with the initial assumption. We can calculate the ordering 



temperature Tb p with the different distributions P(J) derived in section 3.2. We find: 



(i) with the one dimensional model distribution: Tb p = 2(z — l)x£A/(l + 2x£); 

(ii) with the isotropic two dimensional model distribution: Tb p — 2(z — l)-7rx£ 2 A in the dilute 
regime x£ 2 < 1. 

(iii) with the quasi one dimensional model distribution: Tb p — 2(z — l)Trx^ x ^ y A in the dilute regime 

X£ x £y < 1. 

The three limits therefore show a similar behavior Tb p oc (z — l)xA, showing that the short distance 
Wigner repulsion in the spacing distribution Eq. [3| does not affect the ordering properties. Compar- 
ing the Bethe-Peierls ordering temperature to the ordering temperature obtained from the Stoner 
criterion in Ref. [15], we see that z — 1 should be identified with the interchain coupling. The small- z 



regime, relevant to weak interchain correlations, is now discussed in sections 3.4 and 3.5 



3.4 Bethe-Peierls transition with a finite coordination: (i) percolation approxi- 
mation 

We now consider the physics at a finite z = 3. In this regime, the Bethe-Peierls method takes into 
account the inhomogeneities of the magnetization, not included in the Stoner criterion mean field 
solution in Ref. ||15|| . We first consider a "percolation approximation" in which we assume the bonds 
J <C T (J S> T) to be set to zero (infinity) in the effective percolation problem. With z = 3, the 
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Figure 2: The tree structure used in the Bethe-Peierls calculation. The forward branching ratio is 
z — 1. A tree with n + 1 generations is obtained from connecting z — 1 trees with n generations. 



Bethe-Peierls iteration Eq. |4] reads 

X : 



Y tanh (/?J y ) + Ztanh (/?J 2 ) 
1 + YZ tanh (/3J y ) tanh {(5J Z ) ' 



and is approximated into: (i) T <C J y , T <C J^: X ~ Y + Z; (ii) T J y , T <^ J z : X ~ Z; (hi) 
T <C Jy, T S> J 2 : X ~ Y; (iv) T J y , T ^> J z : X ~ 0. The recursion of the average magnetization 
is therefore ((AT)) n+ i ~ 2A((AT)) n , with the percolation parameter 

r+oo 

A = j^ P(J)dJ. (5) 

With the one dimensional distribution, we have A = 1 — (T / A) 2x ^ x , which yields a transition at the 
temperature 

— p(-^), 

exponentially small in l/(x£a;). This behavior is compatible with Ref. |L5| 1 where we have shown 
the absence of a true ordering transition in the percolation approximation of a two dimensional 
anisotropic model, while the model was shown to percolate in a finite size. 

3.5 Bethe-Peierls transition with a finite coordination: (ii) beyond the percola- 
tion approximation 

We now solve the Ising model beyond the percolation approximation. We take into account the 
iteration of small exchanges to lowest order, with the following approximate iteration: (i) T -C J«, 
T < J z : X ~ Y + Z; (ii) T > J y , T < J z : X ~ /3J y Y + Z; (hi) T < J„, T > J z : X ~ Y + /?J Z Z; 
(iv) T Jy, T J z : X ~ [3J y Y + (3J Z Z. The dominant contribution originates from the region 
(iv) of the couplings: «X)) n+ i ~ (2/7>(l - A)((X)) n , with A in Eq. §, and /x = J Q T JP{J)dJ. 
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Gap distribution, Isotropic model 
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Figure 3: Evolution of the gap distribution of the model with isotropic exchanges £ x = £ y = 9, a = 2 
as the temperature is scaled down. The doping concentration is x = 0.01 and the system has a size 
200 x 200. The weight of energy scales ~ T increases as the temperature is decreased. 

With the one dimensional distribution for P(J), we have /x ~ 2x£A and therefore the same critical 
temperature Tb p = 4x£A as in the model with a large connectivity z (with z = 3 in this calculation). 
It is remarkable that the correct treatment of the small exchanges restores a transition temperature 
oc x£A. This shows the relevant role played by energy scales smaller than the temperature. 

The main unsolved question regarding the classical model behavior is to determine whether the 
finite dimensional model has a true thermodynamic transition at a temperature oc x£A. The Bethe- 
Peierls solution orders at a temperature oc x^A because of the strong short range correlations. This 
does not necessarily mean that the finite dimensional model also has a true thermodynamic transition 
at this temperature. Instead, we believe it possible that the classical model has a cross-over to a 
Griffiths physics at a temperature oc x^A and a true thermodynamic transition with a diverging 
correlation length at a temperature T*, which would also be a behavior compatible with a low 
temperature antiferromagnetic susceptibility. At the present stage, we cannot make the distinction 
between these two behaviors. 

4 Quantum isotropic model 

We first consider the artificial situation where the correlation lengths are identical in the two direc- 
tions: £ x = £ y = £ = 9. The tendency to ordering in this isotropic model is overestimated compared 
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(1 = 2, x=0.01, Isotropic model 
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Figure 4: Temperature dependence of the number of effective moments, normalized to the number 
of initial moments xL 2 . We have set £ x = £y = 9, A = 44.7 K, x = 0.01, a = 2. We used square 
systems of dimensions L x L, with L = 50 (O), L = 100 (+), and L = 200 (□). For large system 
sizes, N cS ~ (T/A)xL 2 . 



to the anisotropic model with £ x = 9, fy = 0.1 x £ x . We are lead to consider the class of interactions 



J, 



(-1) 



d x +dy+l 



A exp 



V 



~r- ] + \ ~r- 



6 



(7) 



decaying faster than the interactions in Eq. ^ if a > 1. The cluster RG (see the Appendix) generates 
large energy scales in the parameter range a < ~ 1.2. It turns out that «o < 1 i n the model 
with anisotropic exchanges, and we therefore consider only the regime a > ao — 1-2 in the isotropic 
model. 

The gap distribution is shown on Fig. || for decreasing temperatures. It is visible that the RG 
produces gaps of order of the temperature T unlike in the case of the infinite randomness fixed point 



where the opposite occurs (see Ref. [27] for the one dimensional Heisenberg chain with an infinite 
randomness, random singlet behavior; see Ref. |2q] for the infinite randomness behavior in the two 
dimensional Ising model in a transverse field). As in the Ising model analysis, we calculate the 
susceptibility in two ways: (i) we assume a paramagnetic behavior of the set of effective moments; 
(ii) we incorporate the correlations induced by exchanges A ~ T, in which case an antiferromagnetic 
behavior in the susceptibility is restored. 
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(1=2, x=0.01, Isotropic model 
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Figure 5: Temperature dependence of the the average square effective moment (([S cS ] 2 }} with £ x = 
£ y = 9, A = 44.7 K, x = 0.01, a = 2. We used square systems of dimensions Lx L, with L = 50 (O), 
L = 100 (+), L = 200 (□). The solid line is A/T. The insert shows the low temperature dependence 
of (xL 2 )- 1 ^^} 2 }). 
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Figure 6: Temperature dependence of the Curie constant Tx{T)/x with £ x = £ y = 9, A = 44.7 K, 
and for a = 2. We used square systems of dimension L x L, with L = 50 (O), L = 100 (+), L = 200 
(□). The low temperature regime has been fit to a power law behavior \T jx = 0.39 x T -01 (L = 50, 
solid line), xT/x = 0.375 x T - 135 (L = 100, long dashed line), and X T/x = 0.43 x T - 18 (L = 200, 
short dashed line). The insert shows the low temperature behavior. 
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4.1 Infinite randomness calculation 

We first consider all the exchanges J < T to be set to zero: the set of effective spins is viewed as a 
paramagnet with a susceptibility 

X=TZ2«£ + 1))), (8) 

where iV( eff ) the number of effective spins. We have discarded a prefactor 1/3 in Eq. ||, not relevant to 
the present calculation. The low temperature susceptibility is therefore controlled by two quantities: 
(i) the density of free spins n e ff = ((iV( eff )))/(zi 2 ); and (ii) the magnitude of the effective spin. 

The number of effective moments scales like N e g ~ (T/ A)xL 2 , as it is visible on Fig. ||. The 
squared effective moment shows two regimes: 

(i) High temperature regime: The high temperature average squared effective moment scales like 
(([5 eff ] 2 )) ~ A/T (see Fig. ||). The susceptibility per unit volume is x ~ x/T. 

(ii) Low temperature percolation regime: At low temperature, the squared effective moment scales 
like ((S 2 )) ~ axL 2 , with a some constant (see the insert Fig. ||). The susceptibility per unit 
volume is % ~ (ax 2 L 2 )/ A. In this regime, a cluster has percolated through the finite size 
system. Its magnetization results from summing xL 2 variables Sf = ±1/2, corresponding to 
the two sublattices magnetizations. Therefore, ((S)) ~ VxL 2 , and {{S 2 )) ~ xL 2 . 

The cross-over between these regimes occurs at the temperature scale T co = (A/a)(xL 2 ) -1 , which 
decreases to zero when the system size is increased. Therefore in the thermodynamic limit, only 
the high temperature paramagnetic behavior survives while in a finite size, a low temperature tail is 
present in the susceptibility (see Fig. ||). Now the situation changes when correlations between spins 
coupled by exchanges of order T are included. 

4.2 Finite randomness calculation 

To schematically incorporate the correlations at energy scales of order of the temperature, we con- 
sider as frozen the spins connected by an exchange with a gap between T/2 and T. This freezing 
results in a staggered magnetization because the set of effective moments is unfrustrated (see the 
Appendix). The resulting susceptibility is shown on Fig. [7[ It is visible that is linear in T at 
small T, with therefore a finite susceptibility at zero temperature. This shows qualitatively how an 
antiferromagnetic behavior can be restored because of the correlations at energy scales A ~ T. 

5 Quantum anisotropic model 
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Figure 7: Temperature dependence of the Curie constant xT/x with an isotropic system of size Lx L, 
with L = 100, L = 150. The spins connected by exchanges between T/2 and T have been frozen. 
This results in a finite Curie constant at zero temperature. 

5.1 One dimensional model 

In one dimension, the RG equations of a model in which only AF nearest neighbor exchanges are 
retained can be solved exactly (see Ref. [|l4|). We note x = J/[Max(J)] the exchange normalized 
to the maximal exchange. The distribution of the variable x is P[x) = (//r)x / / r ~ 1 , with T = 
In (A/Max( J)), and f/T = x£/(l + Tx£). The weight on the strongest exchanges x ~ 1 is ~ x£ 
above the cross-over temperature T* = Aexp (— l/x£), and decreases to zero at temperatures below 
T* , where the system has crossed over to the random singlet fixed point. As a test of our program, we 
considered the cluster RG of a one dimensional model, with the exchanges not restricted to nearest 
neighbors (see Eqs. |l|, ||). For any practical temperature above T*, the weight on energy scales of 
order T is found to remain constant as the system is renormalized, with therefore the same behavior 
as in the one dimensional model with AF nearest neighbor exchanges only. 

5.2 Anisotropic model 

We show on Fig. || the evolution of the gap distribution as the system is scaled down, with the 
parameters £ x = 9, £ y = 0.1 x £ x , relevant to CuGeOs. As in the isotropic model, energy scales 
of order T are generated upon renormalizing the system. To qualitatively include the effects of 
correlations at energy scales A ~ T, we consider the spins connected by an exchange with a gap 
between T/2 and T to be frozen, and obtain a low temperature power-law Curie susceptibility, 
X ~ T a , with a = —0.7 if x = 0.01 (see Fig. |9|). The low temperature susceptibility diverges slower 
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Figure 8: Gap distribution of the anisotropic model, with £ x = 9, £ y = 0.1 x £ x , and x = 0.01 and a 
system of size L x x L y , with L x = 640 and L y = 64. The oscillatory behavior at high temperature 
is due to the anisotropy in the couplings. 
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Figure 9: Temperature dependence of the Curie constant xT/x of the anisotropic model, with £ x = 9, 
£ y = 0.1 x £ x , and x = 0.01. The lattice sizes are 320 x 32 and 480 x 48. The solid line is a fit of the 
low temperature behavior to the form %T = 0.14 x T 3 . 



11 



than a Curie law, which is a behavior characteristic of an antiferromagnet. We did not succeed 
to obtain a > as it is the case in doped CuGeC-3. Therefore, we cannot rigorously conclude on 
whether antiferromagnetism is long ranged or associated to a zero temperature transition. A precise 
discussion of this point is an open question, and would require the correlations at energy A ~ T 
to be incorporated beyhond our present treatment. For instance the cluster RG could be used to 
renormalize the high energy physics and the low energy effective Hamitonian could be treated by 
exact diagonalizations. 



6 Conclusions 

We have shown that the physics of the quantum Hamiltonian Eqs. |], ^ was already present at the 
level of the classical Ising model. A Bethe-Peierls treatment of the classical model has been given in 
which a transition at a temperature oc x£A was found. The quantum Hamiltonian has been treated 
in a cluster RG. The model was shown to have a finite randomness behavior. We have shown at a 
qualitative level how a low temperature antiferromagnetic susceptibility can be obtained. 
Two questions are left open: 

(i) The Bethe-Peierls solution orders at a temperature oc x£A. We do not know whether the two 
dimensional model has also a thermodynamic transition at a temperature oc x£A, or whether 
this temperature scale corresponds to a cross-over to a Griffith physics. Both behaviors would be 
a priori compatible with the existence of a maximum in the susceptibility of the antiferromagnet 
at a temperature oc x£A. 

(ii) The quantum model susceptibility shows an antiferromagnetic behavior at low temperature 
due to correlations at energies A ~ T. The isotropic model shows a finite susceptibility at low 
temperature while the quasi one dimensional has a susceptibility diverging slower than a Curie 
law. A precise investigation of the low temperature susceptibility would require a treatment 
going beyond our present analysis, for instance by treating the low energy effective Hamiltonian 
by exact diagonalizations. 

Two other proposals to explain antiferromagnetism in doped in CuGe03 have been made: Fukuyama, 



Tanimoto and Saito and Mostovoy, Khomskii, and J. Knoester |[29| . These proposals are quite 
different from ours, and we have exposed previously why we think they our model is more rele- 
vant [15|. The inclusion of interchain interaction in our model in Ref. and the present work, 



points strongly towards an compatibility with experiments. 
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Figure 10: The first RG transformations in a cluster expansion with a residual spin formation (a) 
and (b), and a singlet formation (c) and (d). The dashed lines represent renormalized exchanges. 

RG calculations as well as the proof that the effective problem remains unfrustrated as the system 
is scaled down. J. Souletie suggested the existence of a similar physics in the classical and quantum 
models. The cluster RG calculations have been performed on the CRAY T3E supercomputer of the 
Centre Grenoblois de Calcul Vectoriel of the Commisariat a l'Energie Atomique. 

A Renormalization equations 

We use a cluster RG to renormalize the quantum Hamiltonian Eqs. ffl, |2[ The method relies on 
a perturbative expansion in the inverse of the largest exchange, and was originally proposed by 
Dasgupta and Ma j3(| in the context of disordered Heisenberg chains. The cluster RG was applied by 
Bhatt and Lee [31| to a model of phosphorus doped silicon. Fisher used the method to solve exactly 



the random singlet fixed point p7| |. The cluster RG was also used to investigate the low energy 
physics of disordered spin chains: the dimerized Heisenberg chain with random exchanges p^]; the 
spin-one chain with random exchanges |33|, p3]; Heisenberg chains with random ferromagnetic and 



antiferromagnetic couplings [35]. Recently, Motrunich et al. [28] shown the existence of an infinite 
randomness fixed point in two dimensions in the Ising model in a transverse field. At such a fixed 
point, inhomogeneities in the disorder grow indefinitely as the system is scaled down, as in the random 
singlet fixed point in one spatial dimension. We now give a short derivation of the RG equations. 
We isolate two spins Si and S2 coupled by an exchange Ji_2- This sets an energy scale given by 
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the gap between the ground state and the first excited multiplet: if J\-2 > is antiferromagnetic, 
the ground state has a spin S = \S± — S2I and the first excited multiplet has S = \Si — S2I + lj with 
a gap Ai_2 = I Ji-2|(|5i — S2I + 1). If J1-2 < is ferromagnetic, the ground state has S = S\ + S2 
and the first excitated multiplet has S = S± + S2 — 1, with a gap Ai_ 2 = | ^1—2 1 ('S'l + Sa)- Among 
all possible pairs of spins, we consider the one with the strongest gap Ai_ 2 - This energy scale is 
identified to the system temperature. If Si and S2 are coupled ferro (J1-2 < 0) , the two spins Si 
and S2 are replaced by an effective spin S = Si + S 2 - If they are coupled antiferro, they are replaced 
by an effective spin S = \Si — S^l- S\ = S2 with an AF coupling Ji_ 2 leads to singlet formation 
while a residual moment is formed otherwise. 

A.l Residual moment formation 

Let us first consider the case where a residual moment is formed corresponding to (a) and (b) on 
Fig. [u]. We specialize a spin S3 among the other spins and denote by Jj_j the exchange between 
spins i and j, with i,j = 1, 3. The coupling Hamiltonian between the spins Si and S2 is -ffi-2 = 
</i-2Si.S2 while the remaining couplings 

H T = J!_3Si.S 3 + J2-3S2.S3 (9) 

are treated in a first order perturbation. This leads to the renormalized coupling Hamiltonian 
Hj = J3S3.S, with the renormalized exchange 

J 3 = Ji_ 3 c(5i,5 2 ,5) + J 2 -sc(S 2 ,S 1 ,S), (10) 

with 

, Q Q q\ SQS + 1)+Si(5i + 1)-S 2 (5 2 + 1) 

c(5l> 52 ' S) = WsTT) 



derived in Ref. [35]. The sublattice on which the residual spin is placed is determined as follows: if 
Si > S2, the residual spin S is placed on the same sublattice as S'l while it is placed on the sublattice 
of S2 if 5*1 < S2. 

A. 2 Singlet formation 

We now consider singlet formation, with S\ = S2 coupled AF. The renormalized couplings are 



obtained in a second order perturbation theory. Generalizing the calculation in Ref. [31, |35| to the 
coupling Hamiltonian @, we find the renormalized exchange 

J3-4 = J3-4 H 7TT ~(^i-3 _ ^2-3X^2-4 — ^1-4), (11) 

3^1-2 



where Si = S2 denote the spins at site 1 and 2. In the ID limit J2-3 = J1-4 = Eq. 11 reproduces 



the result in Ref. [35]| , and the spin- 1/2 limit S\ = 1/2 reproduces the result in Ref. [31]. 
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A. 3 Absence of frustration 



We show that frustration is not generated by the RG procedure. We assume an unfrustrated starting 
Hamiltonian, and show that the different RG operations are compatible with the sublattice structure. 
We distinguish three cases: 

(i) Si and S2 belong to different sublattices and are coupled antiferro. We assume S± > S2 and 
the effective spin S = Si — S2 replaces the spin S\ . The renormalized coupling to another spin 
S3 in Eq. is 

~ S 2 

J3 = J1-3 + (yi-3 — J2-3) g y 

Jl-3 > and J2-3 < leads to J3 > 0. J1-3 < and J2-3 > leads to J3 < 0. The 
renormalized coupling J3 has thus a sign compatible with the sublattice structure. 

(ii) Si and S2 belong to the same sublattices and are coupled ferro. J1-3 and J2-3 have the same 
sign. We have c(Si, S2, S) > and c(Sz,Si,S) > 0. The renormalized coupling J3 has the 
same sign as J1-3 and J2— 3, compatible with the sublattice structure. 

(iii) Si = S2 are coupled antiferro and a singlet is formed. If S3 and S4 are coupled ferro and in 



the same sublattice, J3-4 < in Eq. |11|. If S3 and S4 are coupled antiferro and in the opposite 
sublattice, J3-4 > in Eq. 11. The singlet formation is thus compatible with the sublattice 
structure. 
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